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In this paper, we will describe the results of SUSY parameter space searches including minimal supergravity, 
non-universal supergravity and minimal supersymmetry and the implications on the indirect detection of neu- 
tralino dark matter. We give special attention to the effects of detector thresholds, solar absorption of neutrinos 
and hadronization of neutralino annihilation products. These effects are known to be important in calculating 
accurate event rates [El. 

We chose also to focus on models which predict a heavy lightest neutralino (several hundred GeV to several 
TeV). These models have been selected for several reasons including their inaccessibility in future collider searches. 



1. INTRODUCTION 

There has been a great deal of evidence mount- 
ing recently for the existence of dark matter con- 
sisting primarily of the the lightest supersymmet- 
ric particle (LSP) [ |j. In most models, the LSP 
is the neutralino, a superposition of the super- 
symmetric partners of the photon, Z boson and 
neutral Higgs bosons. By virtue of i?-parity con- 
servation, the LSP is totally stable and would 
freeze out during early stages of the universe. In 
many models, the predicted relic density of neu- 
tralinos is in the cosmologically interesting region 
(.05 < Vl x h 2 < .3) which is further evidence for a 
significant neutralino dark matter component. 

In this paper, we consider a large area of min- 
imal supergravity (mSUGRA), non-universal su- 
pergravity and minimal supersymmetry (MSSM ) 
parameter spaces. We determine which regions 
predict an acceptable cold dark matter particle 
candidate and calculate the upward going muon 
fluxes in neutrino telescopes from neutralino an- 
nihilation in the sun and compare these results 
with direct detection rates for these models. We 
have also calculated the rates from annihilation in 
the Earth, but do not discuss these results here 
as they are generally negligible compared to the 
solar rates. 

The results presented in this paper are de- 
scribed in more detail in Reference [ 01 . 



2. BACK-OF-THE-ENVELOPE CALCU- 
LATIONS 

In this section, we roughly calculate the event 
rates for indirect and direct WIMP detection ex- 
periments in a model independent fashion. No 
references to any specific supersymmetry or other 
model are required for the following calculations. 
These calculations closely follow Reference [ Q . 

2.1. Assumptions 

Before starting the calculations, some basic as- 
sumptions need to be stated. First, we assume 
that the WIMP in question is the major con- 
stituent of the measured dark matter halo den- 
sity: 



(1) 



We will also approximate, using dimensional 
analysis, the cross section for WIMP-Nuclcon in- 
teractions: 
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Thirdly, we will assume that WIMPs annihi- 
late to neutrinos ~10 percent of the time. This 
is simply a reflection of the branching ratios of 
heavy quarks or W boson pairs to neutrinos. 
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2.2. Indirect Detection 

Indirect detection experiments are based on the 
idea that heavy particles, such as WIMPs will be- 
come gravitationally bound in objects such as the 
Sun and Earth. Over a time scale generally less 
than the age of the solar system, these bodies can 
obtain WIMP densities which will produce an ob- 
servable flux of annihilation products (generally 
neutrinos) for a detector on Earth. 

With the assumptions stated and out of the 
way, we will begin the calculation of indirect event 
rates from solar WIMP annihilation. First, the 
solar gravitational capture cross section is given 
by: 

a sun = f(1.2 x 10 57 ) auA (3) 

Where f is the focusing factor and is equal to 
the ratio of kinetic to potential energy of the 
WIMP near the Sun and is generally of the order 
of 10. The flux of WIMP annihilation products 
from the Sun in units of number of WIMPs is: 

</> SU n = (/) x a sun /4ird 2 (4) 

Where d is the Earth-Sun distance. Then, us- 
ing the assumption of ten percent branching ra- 
tios of WIMPs to neutrinos, the resulting neu- 
trino flux at Earth is given by: 

3 x 10 -5 

4> v = .10 x <f> snn = — — cm^s" 1 (5) 

m x l GreV 

The probability to detect a neutrino in an ice- 
based neutrino telescope is proportional to the 
neutrino interaction cross section and the muon 
range within the detector. This, along with the 
assumption that the neutrino energy is half (or 
one third in the case of heavy quark decay) of the 
WIMP mass, gives us the probability of detecting 



a neutrino as: 

P = N A a u R^ (6) 

a v = 0.5 x 10- 38 (£J„/ GeV) cm 2 (7) 

i? A1 =muon range= 500 cm (E M / GeV) (8) 

-> P = 2 x 10~ 13 (m x / GeV) 2 (9) 



Thus, the indirect event rate is simply: 

^indirect - <f>»P = 1.8 (m x /GeV) y^km" 2 (10) 

This calculation is very approximate in na- 
ture and is only intended to provide a rough 
idea of the factors involved in predicting such an 
event rate. Some of the neglected factors include 
the solar absorption of neutrinos, the hadroniza- 
tion of quarks, neutrinos from other annihila- 
tion channels, detector threshold effects, reso- 
nances, model-dependent characteristics of the 
WIMP, coannihilations and coherent nuclear en- 
hancements. Later in this paper, most of these 
factors will be considered and included in the re- 
sults given. 

2.3. Direct Detection 

For the purpose of comparison, we will now pro- 
ceed to calculate, in a similar manner, the event 
rates for direct detection experiments. Direct ex- 
periments measure the energy transferred from 
halo WIMPs directly to the experimental appa- 
ratus. The event rate for such an experiment is 
given by: 

^Direct = — 4>X a DA = TTT-TT kg _1 yr _1 (H) 

Tji 7v m x / GreV 

This result provides a rough estimate of such an 
event rate. It neglects, among other factors, the 
model-dependent characteristics of the WIMP, 
coherent nuclear enhancements and the proper- 
ties and composition of the detector. 

Comparing these approximate results for in- 
direct and direct rates, we can gauge the rel- 
ative effectiveness of such experiments. For a 
10,000 square meter indirect detector and a one 
kilogram direct detector, direct detectors observe 
more events only for the case of an extremely light 
WIMP (less than 10 GeV) which has previously 
been ruled out [ |) . This result is not general due 
to the factors neglected in the above calculations. 
A much more detailed comparison will be given 
later in this paper. 

In addition to signal comparison, the back- 
ground for indirect and direct detection experi- 
ments should be considered. Indirect detection 
background is from atmospheric neutrinos and is 
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Figure 1. Effects of solar absorption, hadronization and threshold on indirect event rates. The models 
considered in the figure are the MSSM with = 400 GeV, tan (3 = 35 and reflect typical behavior for 
most models with a similar mass. Heavier models are effected much less by such factors. 



determined by their known energy spectrum. The 
number of events can vary from a few thousand 
background events per year per solar pixel for 
neutrinos of energy around 10 GeV to less than 
one background event per year per solar pixel 
for neutrinos of energy near 1 TeV. Direct de- 
tector backgrounds are primarily energy indepen- 
dent and are generally around 300 events per year 
per kilogram of detector material. Therefore, at 
high energies, indirect detectors have background 
advantages while at low energies, direct detectors 
have fewer background events. 

3. Effects of Hadronization, Solar Absorp- 
tion and Detector Threshold on Indirect 
Detection Rates 

To remove as much of the uncertainty in the 
above calculation as possible, we have taken steps 
to include many additional factors in our pa- 
rameter space searches. To do so, we first con- 



sider the dominant annihilation modes for a neu- 
tralino WIMP. For a gaugino-like neutralino, the 
dominant annihilation modes are bb or rf for 
m x < m w ± and bb for m x > m w ± . For a 
higgsino-like neutralino, the dominant modes are 
ZZ or W + W~ for m x < m t and ti for m x > m t 
[ 0). Any quark products, and the resulting neu- 
trino products, lose a fraction of their energy to 
hadronization as determined by the quark flavor 




Furthermore, neutrinos produced in the Sun 
may be absorbed before they escape the solar 
medium. We have used the Ritz and Seckel pa- 
rameterization up to energies around 1 TeV for 
all neutrino producing annihilation modes [ 
Above the TeV scale, corrections were needed and 
were done in a manner similar to work by Edsjo 




Indirect detectors generally have neutrino en- 
ergy thresholds of a few tens of GeV. We have cho- 
sen 25 GeV for the results in this paper. We have 
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Figure 2. Indirect detection rate of neutrinos from solar WIMP- WIMP annihilation in the MSSM. Regions 
represent more than 100, more than 10 and less than 10 events per year per square kilometer effective 
area (darkest to lightest). 



also made the calculations for a threshold of 100 
GeV but in the interest of space have neglected 
these results. For a more detailed description of 
these effects, see References [ ||, ED]. 

The total effect of hadronization, solar absorp- 
tion and threshold can vary with energy from the 
total destruction of the neutrino signal for masses 
around 100 GeV to merely a reduction by a fac- 
tor of two or less in neutrino signal at TeV scale 
masses. Figure 1 shows the results of these effects 
for a typical SUSY model with a light neutralino 
of about 180 GeV. The effects are far less signifi- 
cant for heavier models. 

The calculations presented in this discussion 
neglect the effects of coannihilations. These ef- 
fects are presented in Reference [ |) . 

4. Searches of Supersymmetric Parameter 
Space 

Beginning with the most simple set of SUSY 
models, the MSSM, we attempt to identify the 



model-dependent properties of neutralino dark 
matter detection. To simplify the parameter 
space, we set all squark and slepton masses to 300 
GeV or 1.5 times the neutralino mass (whichever 
is larger) [ EH]]. We have also made calculations 
for the choice of 1500 GeV or 1.5 times the neu- 
tralino mass but in the interest of space do not 
display them here. We then search by varying 
fi, M 2 , tan/? and the couplings. We set Mi rel- 
ative to M 2 using the GUT scale relations and 
we hold M H2 ~ 120 GeV (within 15 GeV). We 
consider only those choices of parameters which 
yield a neutralino with a cosmologically interest- 
ing relic density (.05 < il x h 2 < .3) which is the 
lightest supersymmetric particle (LSP) . Figures 2 
and 3 show the indirect detection rate and ratio of 
indirect to direct detection rates for the MSSM. 
Note that when the lightest neutralino is largely 
gaugino, which is most often the case, the GUT 
relations cause m x to scale with M 2 . 

Extending the MSSM to minimal supergravity 
(mSUGRA), we simplify our model parameters 
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Figure 3. Ratio of solar indirect detection rate to direct detection of galactic halo WIMPs in the MSSM. 
Regions represent more than 100, more than 10 and less than 10 in units of events per year per square 
kilometer effective area and events per year per kilogram of detector material (darkest to lightest). 



to the universal scalar mass at the GUT scale 
mo, the universal gaugino mass the GUT scale 
mi/2; the ratio of vacuum expectation values of 
Higgs scalars tan (3, the universal trilinear cou- 
pling at the GUT scale Aq and the sign of the 
superpotential term in the Lagrangian fi. Re- 
sults shown here are for models with ^4o=0 and 
(j, > 0. We have calculated results for other mod- 
els and found generally no significant differences. 
Changes in A mostly effect masses of third gen- 
eration superpartners and, therefore, have little 
effect on the relic density or relevant cross sec- 
tions [ Figures 4 and 5 show the indirect 
detection rate and ratio of indirect to direct de- 
tection rates in mSUGRA. 

Extending supergravity to the cases of non- 
universality, we introduce two new parameters S\ 
and S 2 which parameterize the divergence of the 
Higgs masses from the universal scalar mass at 
the GUT scale: H x = {l+5 1 )m ,H 2 = {l + S 2 )m Q 
[ [l§. The case of S% = S 2 = reduces to 
mSUGRA. In both minimal and non-universal 



SUGRA, we consider only models which allow for 
acceptable electroweak symmetry breaking and 
are tachyon-free in addition to the requirements 
imposed above in the MSSM models. Figures 6 
through 9 show the indirect detection rate and 
ratio of indirect to direct detection rates for the 
non-universal SUGRA. Note that the mass of the 
lightest neutralino scales with mw2 and is roughly 
40% of m X j 2 for a largely gaugino lightest neu- 
tralino, which is the case for most mSUGRA and 
non-universal SUGRA models [ [b| [b| 

5. Discussion 

Figures 2 and 3 show the indirect detection 
rates and the ratio of indirect to direct detection 
rates for the MSSM. The regions represent more 
than 100, more than 10 and less than 10 (in units 
shown on the figures). Empty regions represent 
models with non-cosmologically interesting neu- 
tralino relic densities. 

The rates are clearly more favorable for tan (3 — 
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Figure 4. Indirect detection rate of neutrinos from solar WIMP- WIMP annihilation in mSUGRA. Regions 
represent more than 10, more than 1 and less than 1 in units of events per year per square kilometer 
effective area (darkest to lightest). Theoretically disfavored regions are discussed in the text. 



35 than for tan j3 = 3 models. It is difficult to say 
how many events are necessary for a positive dark 
matter detection. A fair estimate of the number 
needed for a positive result is probably around 
10 events per year for a model with a lightest 
neutralino with mass around a few hundred GeV. 
With this in mind, much of MSSM parameter 
space should be experimentally accessible to kilo- 
meter scale neutrino telescopes. From Figure 3, 
it is clear that indirect detection is more favor- 
able than direct detection in the vast majority of 
MSSM models. This advantage is even greater, in 
all but the lightest models, if the effects of back- 
ground are included. 

The searches of mSUGRA and non-universal 
SUGRA space are somewhat less favorable than 
the MSSM but still show promise. The regions 
in Figures 4 through 9 represent more than 10, 
more than 1 and less than 1 (in units shown on 
the figures). Empty regions are either disfavored 
by their neutralino relic density or by other theo- 
retical constraints such as acceptable electroweak 



symmetry breaking, the existance of tachyons. 

Much of the SUGRA parameter space produces 
a few events per year for a kilometer scale neu- 
trino telescope. This would likely require many 
years of observation to observe a positive result. 
Large portions of the tan (3 — 50 region would not 
be observable even with many years of observa- 
tion. 

6. Conclusions 

Several important conclusions can be drawn 
from the results presented in this paper. First, 
in most SUSY models, indirect detection is fa- 
vored over direct detection methods, especially 
in heavier neutralino models. Second, the effects 
of hadronization, solar absorption and detector 
threshold are significant and should not be ne- 
glected. Third, a significant portion of super- 
symmetric parameter space, especially within the 
MSSM should be accessible to present or next 
generation indirect detection experiments. More 
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Figure 5. Ratio of solar indirect detection rate to direct detection of galactic halo WIMPs in mSUGRA. 
Regions represent more than 10, more than 1 and less than 1 in units of events per year per square 
kilometer effective area and events per year per kilogram of detector material (darkest to lightest). 
Theoretically disfavored regions are discussed in the text. 



importantly, however, is that the areas of SUSY 
parameter space which are generally most accessi- 
ble to indirect detection are the same areas which 
will be very difficult to probe in collider experi- 
ments [ [l(| [l?], [l]|. With this in mind, if na- 
ture has selected a heavy scale for supersymme- 
try, indirect dark matter experiments may make 
the first observations of supersymmetry. 
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Figure 6. Indirect detection rate of neutrinos from solar WIMP- WIMP annihilation in the non-universal 
SUGRA case of Si = —.5. Regions represent more than 10, more than 1 and less than 1 in units of events 
per year per square kilometer effective area (darkest to lightest). Theoretically disfavored regions are 
discussed in the text. 
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Figure 7. Indirect detection rate of neutrinos from solar WIMP- WIMP annihilation in the non-universal 
SUGRA case of #2 = -5. Regions represent more than 10, more than 1 and less than 1 in units of events 
per year per square kilometer effective area and events per year per kilogram of detector material (darkest 
to lightest). Theoretically disfavored regions are discussed in the text. 
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Figure 8. Ratio of solar indirect detection rate to direct detection of galactic halo WIMPs in the non- 
universal SUGRA case of 5i = —.5. Regions represent more than 10, more than 1 and less than 1 in 
units of events per year per square kilometer effective area (darkest to lightest). Theoretically disfavored 
regions are discussed in the text. 




Figure 9. Ratio of solar indirect detection rate to direct detection of galactic halo WIMPs in the non- 
universal SUGRA case of 82 = -5. Regions represent more than 10, more than 1 and less than 1 in units of 
events per year per square kilometer effective area and events per year per kilogram of detector material 
(darkest to lightest). Theoretically disfavored regions are discussed in the text. 



